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We report measurements of the inclusive electron momentum spectra in decays of charged and 
neutral B mesons, and of the ratio of semileptonic branching fractions B(B + —* Xev) and B(B° — > 
Xev). These were performed on a sample of 231 million BB events recorded with the BABAR 
detector at the T(45) resonance. Events are selected by fully reconstructing a hadronic decay of 
one B meson and identifying an electron among the decay products of the recoiling B meson. We 
obtain B(B+ -> Xev)/B(B° -» Xev) = 1.084 ± 0.041 (stat) ± 0.025 (syst) . 

PACS numbers: 12.15.Hh, 11.30.Er, 13.25.Hw 



The hadronic decay widths of B + and B° mesons dif- 
fer because of mechanisms that depend on the flavor 
of the spectator quark, such as interactions involving 
the spectator quark or final state particles. This leads 
to different lifetimes t b + and t b o of charged and neu- 
tral B mesons. We do not expect different semilep- 
tonic decay widths, since semileptonic decays do not in- 
volve the spectator quark. This means that the ratio 
R +/0 = B(B+ -> Xev)/B(B° -> Xev) should agree with 
Tg+ /Tgo , which can be checked experimentally. 

At the T(45 l ) resonance, measurements of the inclusive 
semileptonic branching fractions of B + and B° mesons 
are less precise than for an admixture of b hadrons. The 
reason is mainly a limitation of statistics from the small 
efficiency of the event tag needed to separate B + B~ from 
B°B° events. In this paper, we use fully reconstructed 
hadronic B decays for this separation. Combined with 
the high statistics of the B factories, this approach al- 
lows for a precision measurement of R+/o, as already 
demonstrated by the Belle collaboration, measuring R + / 
with 5% uncertainty 0. By tagging B°B° events with 
partially reconstructed B° — > D*~ Iv decays, the CLEO 
collaboration achieved a 14% uncertainty on R+/o 
High-momentum electron tags have been used in similar 
analyses for the determination of B(B — > Xev) and the 
electron momentum spectrum without separation of B° 
and B+ decays 0,0]. 

The measurements presented here are based on data 
collected by the BABAR detector 5] at the PEP-II asym- 
metric e + e~ storage rings and correspond to an inte- 
grated luminosity of 209 fb _1 (231 million BB events) 
on the T(45) resonance. For background and effi- 
ciency corrections that cannot be measured directly from 
data, we use a full simulation of the detector based on 
GEANT4 0- The equivalent luminosity of the sim- 
ulated event sample amounts to about 980 fb _1 for 
T(4S f ) — > BB events and 300 fb _1 for processes consist- 
ing of non-resonant e + e~ — + qq (q = u, d, s, c) production 



("continuum events"). 

In events with a fully reconstructed hadronic B decay 
(Stag), we identify electrons among the remaining tracks. 
To avoid large backgrounds at lower momenta, we re- 
quire p e > 0.6 GeV/c, where p e is the electron momentum 
measured in the center-of-mass frame. Depending on the 
electron charge q e relative to the charge qb of the bottom 
quark in the -B tag candidate, each electron is assigned to 
either the right-sign (q e — —iqb) or to the wrong-sign 
sample (q e — 3qb). In events without B°B°-mixing and 
a correctly reconstructed Bt ag , primary electrons from 
semileptonic decays of the signal B are the dominant 
source for the right-sign sample, while electrons from 
B — > DX,D — > e~v e Y cascades populate the wrong- 
sign sample. We use the criteria in Ref. [3j for track 
selection and electron identification, and apply the same 
procedures for efficiency and background corrections of 
the right- and wrong-sign samples. In this analysis, we 
additionally have to correct for mis-reconstructed B tag 
candidates. 

Non-BB events are suppressed by requiring the ratio 
of the second to the zcroth Fox- Wolfram moments 
to be less than 0.5. To keep backgrounds from mis- 
reconstructed B tag candidates at a low level, we recon- 
struct hadronic B decays in very pure modes only. To 
cancel systematic errors related to the B tag reconstruc- 
tion, we select similar ("twin") modes for B° and B + 
decays 
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Here tt and K denote charged pions and kaons. The 
invariant mass of D° candidates is required to be within 
15 MeV/c 2 of the nominal D° mass |9( for the decay 
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D a -> Kir and 25 MeV/c 2 for D° -> i^Trvr decays. D" 
candidates are accepted if the invariant mass is within 20 
McV/c 2 of the nominal D~ mass. D candidates with mo- 
menta above 2.5 GeV/c (measured in the center-of-mass 
frame) are rejected since they indicate non-BB events. 
D* candidates are built from pairs of D° candidates and 
charged (neutral) pions where the invariant mass differ- 
ence \Mjjo n (o) — Mjj \ is within 2 MeV/c 2 of the nominal 
mass difference. In tag categories (III) and (V) we re- 
quire the invariant masses and M wn o to be less 
than 1.5 GeV/c 2 . For further background reduction, we 
reject candidates where a kinematic fit with geometric 
constraints on the B and D vertices and mass constraints 
on the charmed mesons yields a \ 2 value with a proba- 
bility of less than 0.5%. 

The kinematic consistency of the Stag candidates is 
checked with two variables, the beam-energy substituted 
mass tubs = 

(s/4 - p%) 1/2 and the energy difference 
AE = Eb — y/s/2. Here y/s refers to the total center- 
of-mass energy, and Eb and ps denote the energy and 
momentum of the Btag candidate, all quantities being 
measured in the center-of-mass frame. For categories 
(I) - (III), we require \AE\ < 50MeV, while the pres- 
ence of an additional 7r° in (IV) and (V) leads to asym- 
metric distributions in AE, motivating lower limits of 
AE > -75MeV for (IV) and AE > -100 MeV for (V). 
If for a given mode more than one -B tag candidate sat- 
isfies these criteria, the one with the smallest \AE\ is 
selected. Figure ^ shows the tues distributions of B tElg 
candidates satisfying these selection criteria. Candidates 
with 5.27 < niES < 5.29 GeV/c 2 are included in the B tag 
sample. In w 1% of all events, we find multiple -B tag can- 
didates in different decay modes. Here we use all of them, 
correcting for the background Stag candidates later. 

The -Btag sample can be divided into 4 components: 
signal, combinatoric background, D*~ «-> D*° cross feed 
and continuum background. Correctly reconstructed B 
decays are called signal Btag candidates, while -Btag can- 
didates that contain tracks from the decay of the other 
B contribute to the combinatoric -B tag background. A 
special case of combinatoric background, called D*~ <-> 
D*° cross feed, contains cross feeds between twin modes 
of channels (II)-(V) due to mis-reconstruction of a D*~ 
as a D*° or vice versa. Due to the low energy of the com- 
binatoric pion, the uies distribution of this background 
is similar to the signal and will be treated separately from 
the other combinatoric -B tag background. The fourth 
component consists of -B tag candidates arising from con- 
tinuum events and is called continuum -B tag background. 
Since the ratio of signal to background -B tag candidates 
depends on the multiplicity of the event and thus on the 
presence of a semileptonic decay, a precise determina- 
tion of the number of signal -B tag candidates is crucial 
to avoid biases in the branching fraction measurement. 



Monte Carlo (MC) studies using generator information 
indicate that once the -B tag , right- and wrong-sign sam- 
ples have been corrected for -B tag background, the biases 
on the branching fraction measurements are below the 
statistical sensitivity given by the size of the MC sample, 
i.e. less than 0.5%. 
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FIG. 1: Fits of Eq.Qto distributions of the energy substituted 
mass for (a) neutral and (b) charged -B tag candidates. The 
dotted and dashed curves indicate the fitted contributions of 
continuum and combinatoric -Btag candidates. The grey his- 
togram displays the contribution of D*~ <-» D*° background. 



The contributions of combinatoric and continuum -B tag 
background to the -B tag sample are extrapolated from the 
niEs sideband region, 5.2 < tues < 5.25GeV/c 2 . This 
requires a model of the background ities distributions 
over the full range, 5.2 < m E s < 5.29GeV/c 2 , which is 
obtained by fitting a linear combination of three func- 
tions describing the shapes of tties distributions of sig- 
nal, combinatoric and continuum -B tag candidates to the 
observed tues distributions. 

The shape of the combinatoric -Bt ag background 
fbb( mES ) i s taken from the MC simulation. For the con- 
tinuum background, we use the following function |Toj : 



fqq( m ) 



where m = rriEs/ m ES X ano - m Es yi ^ s the endpoint of the 
ities distribution. 

For a given B decay mode, the signal tues distribu- 
tion is commonly described by a gaussian and a power 
law • Since the -B tag signal consists of many individual 
decay modes, a single function of that type fails to de- 
scribe our ities distribution. We have found that a more 
general ansatz using a gaussian shape f g {x) — e~ x I 2 and 
a function with a similar shape near x = 0, but behaving 
like e~ x for x — ► ±oo, ft(x) = e~ x /(l + e~ x ) 2 , yields a 
good description of our signal mES shape: 



C-2 



/»ig(A) = 



(C 3 -A)» 

1 <?L 



if A < 



ifa< A<0 , (1) 
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with A = ities — mES and wbs being the maximum of 
the ities distribution. C\, C2 and C3 are functions of 
the parameters Mes, r, ax, a 2 , <Jl, ol and n to ensure 
that / s ig is continuous and differentiable at A = and 
A = a. This function, similar to the one featured in [ill] , 
describes the tails caused by the asymmetric energy reso- 
lution of neutral pions by a power law of order — n and a 
junction a < where it turns into a gaussian-like shape. 
Fixing a and n to the values obtained from a fit to MC- 
simulated vtles distributions of signal -B tag candidates, 
we fit a linear combination of f q -q, f b ^ and / s j g to the tues 
distributions observed in data, leaving all other param- 
eters and normalizations free in the fit (Fig.^l. Due to 
their similar ities distributions, this method cannot dis- 
tinguish between signal B tag candidates and D*~ <-> D*° 
cross feed. This background contribution is estimated 
from the MC simulation to be 0.5% (2.6%) relative to 
the signal for the neutral (charged) B tElg sample. 

To validate this extraction method, we perform the 
same analysis on our Monte Carlo sample and find that it 
reproduces the original number of signal Stag candidates. 
Uncertainties related to the MC simulation of the combi- 
natoric B tag background are evaluated by decomposing 
this background into the true underlying individual ex- 
clusive decay modes, and varying their contributions by 
the uncertainties of their branching fractions if they are 
reported in Q, or ±100% otherwise. This leads to an 
uncertainty of 1.3% on the number of B° and B + tags. 
Due to the different compositions of the combinatoric B° 
and B + backgrounds, these errors are uncorrelated. In 
contrast, systematic errors related to the description of 
the signal shape are correlated since we use similar de- 
cay modes. Here we assess the uncertainties related to 
the modeling of the shape for tubs < ttles by repeat- 
ing the fit with a set to -00, allowing an exponential 
function only instead of a power law to describe the tail 
caused by the ir° energy resolution. This leads to rela- 
tive uncertainties of 2.1% (2.4%) on the number of B° 
(B + ) tags. The yields of events in which B tag candi- 
dates have been found for both "twins" of decay chan- 
nels (II)- (V) differ by 20% in data and MC, motivating 
a relative uncertainty of 20% on the D*~ <-* D*° cross- 
feed. This adds another systematic uncertainty of 0.5% 
to the number of charged -B tag candidates. The fi- 
nal numbers of neutral and charged signal -B tag candi- 
dates are N B o = 45420 ± 420 (stat) ± 591 (u) ± 949 (c) and 
N B + = 41948 ± 463 (stat) ± 596 (u) ± 1020 (c) , where u and 
c denote uncorrelated and correlated systematic uncer- 
tainties respectively. 

The requirement of an identified electron leads to sig- 
nificantly lower Stag backgrounds, as shown in Fig. [21 
for the right-sign sample. For high electron momenta 
(p e > lGeV/c), the purities are 96% (98%) for the 
right-sign (wrong-sign) samples, with combinatoric -Btag 



candidates being the dominant background, while for 
decreasing electron momenta, the purities decrease to 
90% because of an increasing amount of continuum- 
background. As for the full -Btag sample, we estimate 
these backgrounds from the titles sideband region. The 
background estimates are performed separately for each 
sample as functions of p e . Due to low statistics, we do 
not determine the extrapolation factor from a fit, but 
use the MC predictions instead. The systematic errors 
due to the shape of the combinatoric background and 
D*~ <-> D*° cross feed are evaluated in the same way as 
for the Btag sample, and the uncertainty in the contin- 
uum contribution is taken to be 20%. 
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FIG. 2: ities distributions for (a) neutral and (b) charged 
Stag candidates in events with a right-sign electron. 



Figure shows the momentum spectra of right- and 
wrong-sign electrons in events with a charged -Bt ag candi- 
date, together with the estimated -B tag background. This 
figure also displays the background contributions of elec- 
trons from photon conversions, ir° — > je + e~ Dalitz de- 
cays and misidentified hadrons. These backgrounds are 
identified and corrected for as in 0, . Corrections for 
electron identification efficiency and the evaluation of its 
systematic uncertainty are also performed as in 0,0. 

Background contributions from decays of charmed 
mesons produced in b — > ccs decays or decays of r lep- 
tons are estimated from the MC simulation, using the 
ISWG2 model 12] to describe semileptonic D and D s 
meson decays. Assuming T(D S — > Xev) — Y(D — ► Xev), 
we obtain B(D S -> Xev) = (8.05 ± 0.66)%. Inclu- 
sive D s production has been measured in [l3j ser 
rately for neutral and charged B decays, and with 
we obtain B{B° -> L>+ -> e+) = (0.67 ± 0.17)% and 
B(B+ -> D+ -> e+) = (0.88 ± 0.18)%. Combining 
the measurements of inclusive D° and D + production 
from [l3j with the inclusive D 0,+ — > e branching fractions 



yields B(B° 
,0 



D +,o 



(0.83 ±0.25)% and 



from [_. 

B(B+ ^ £>+<° ^ e+) = (1.33 ± 0.20)%. Since there are 
no branching fraction measurements for B — * t decays 
that distinguish between neutral and charged B decays, 
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FIG. 3: Total measured spectrum (points) and estimated 
backgrounds (histograms) for electron candidates in events 
with a charged Stag candidate, for (a) the right-sign sample, 
and (b) the wrong-sign sample. 



we assume T(B° — > Xtv) = T(B + — > Xtv) and com- 
bine the average value from || with the B-meson life- 
times from direct measurements 0. Including r leptons 
that originate from B — > D s — > r cascades, we arrive at 
B(B° -> r -> e+) =(0.53±0.06)% and £(B+ -> r -> e+) 
= (0.60 ±0.06)%. Since the branching fractions of B de- 
cays to J/tp and ip(2S) mesons are small and well mea- 
sured, we use the MC simulation to correct for back- 
ground electrons from J/ip — > e + e~ and iJj(2S) — > e + e~ 
decays, using B(B -> J/ijj -> e+e") = (6.49±0.22) x 10" 4 
and B(B -> ip(2S) -> e+e") = (0.23 ± 0.02) x 10" 4 0. 

After all corrections listed in Table [I] have been ap- 
plied, the inclusive momentum spectrum of electrons 
from semileptonic decays of B + mesons dN B +^>xev/dp 
is given by the right-sign sample in £?~-tagged events. 
Because of B°B° oscillations, electrons from B° — > Xev 
decays and B° — > DX,D — > e~^ e F cascades contribute 
to both momentum spectra dN r ^ g /dp and dN^/dp of 
right- and wrong-sign samples in £?° - tagged events, 

V 1 Am J r 



dp 



dp 



diV^ _ diV B o^ Xe „ 



dp 



dp 



X™ + 



dp 
dp 



(1 - Xm) , 



with Xm = (0T86 ± 0.004) being the mix- 
ing parameter. We use these equations to determine 
the primary electron spectrum dN B —>Xev / 1 dp of neu- 
tral B decays, which is shown in Fig. 0] together with 
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FIG. 4: (a) Normalized momentum spectra of primary elec- 
trons after all efficiency corrections and (b) their ratio r + / 
= N B o/N B + (dN B +_ Xe Jdp) I (dN B0 _ Xe Jdp). 



dN B +^,xev/dp after normalizations to the respective 
number of tags. 

We integrate these spectra between p min = 0.6 GeV/c 
and 2.5 GeV/c and apply corrections for geometrical 
acceptance (e gcom = 85%) and the small loss of elec- 
trons due to bremsstrahlung in the detector material 
(fbrem = 97.4±0.1%) to obtain the partial branching frac- 
tions B(B° -> Xev{i)) = B(B° -> Xeu(rf),p e > p min ) 
for decays with any number of photons in the final state: 

B(B° Xev(i)) = (9.66 ± 0.27 (stat) ± 0.32 (syst) )%, 
B(B+ Xev(j)) = (10.39 ± 0.26 (stat) ± 0.37 (syst) )%, 
B(B -> Xev(i)) = (10.03 ± 0.19 (stat) ± 0.32 (syst) )%. 

Table |n] lists the contributions to the systematic errors. 
These results are in agreement with y, |4|. For the 
ratio of branching fractions, R + / Q (p m i n ) = B(B + — > 
Xev(~f),p e > p mhl )/B(B° — > Xev(-~f),p e > p mi „), the re- 
sult is i? +/0 (0.6GeV/c) = 1.076 ±0.040 (sta t)±0.029 (syst) . 
For higher values of p m i n , the statistical error increases, 
while the systematic error decreases. At p m in = 1 GeV/c, 
the combined statistical and systematic error is minimal, 
leading to our final result 

i? +/o (1.0GeV/c) = 1.084 ± 0.041 (stat) ±0.025 (syst) . 

In summary, we have used electrons in T(4S*) decays 
tagged by a fully reconstructed hadronic B decay to 
measure the inclusive semileptonic branching fractions 
of B° and B + mesons. The ratio of branching frac- 
tions, A f/0 (1.0GeV/c) = 1.084 ± 0.048, is consistent 
with the ratio of B lifetimes from direct measurements, 
t b +/t b o = 1.086 ±0.017 [9j. From this we conclude that 
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TABLE I: Electron yields for the four samples and corrections with statistical and systematic errors. 





B° tags, 


B° tags, 


B + tags, 


B + tags, 




right-sign 


wrong-sign 


right-sign 


wrong-sign 


5.27 < m BS (B tag ) < 5.29 GeV/c 2 3461 ± 59 


1943 ± 44 


4074 ± 64 


1070 ± 33 


B ta g background 


198 ± 16 ± 40 


135 ±13 ±27 


320 ± 24 ± 64 


114 ±12 ±23 


7 — ¥ e + e~ 


55 ± 14 ± 8 


87 ± 17 ± 12 


66 ± 14 ± 10 


83 ± 16 ± 11 


7T° — > je + e~ 


31 ±14 ±7 


25 ± 12 ± 5 


36 ±14 ±7 


47 ± 16 ±9 


fake e 


29 ± 1 ± 8 


21 ±1 ±4 


37 ± 1 ± 12 


16 ± ± 2 


Yield before and 


3149 ± 64 ± 42 


1674 ± 51 ± 30 


3616 ± 71 ± 66 


810 ±41 ±27 


after e efficiency correction 


3443 ± 70 ± 71 


1842 ± 56 ± 50 


3947 ± 78 ± 96 


898 ± 46 ± 41 


B (A, -»)t e 


92 ± 9 ± 8 


20 ± 4 ± 2 


109 ± 10 ± 9 





S -f D s -» e 


65 ± 9 ± 16 


13 ± 4 ± 3 


96 ± 11 ± 20 





B -> D -> e 


61±8±25 


12 ± 4 ± 5 


96 ± 11 ± 15 





B J/ip,ip(2S) -> e 


22 ± 5 ± 1 


23 ± 5 ± 1 


17±4± 1 


18±4±1 


D »- ^ ^.o crosg feed 


9±3±5 


4±2±2 


44 ± 7 ± 22 


29 ± 5 ± 15 


Net e yield 


3195 ± 72 ± 82 


1769 ±57 ±51 


3585 ± 81 ± 106 


850 ± 47 ± 45 



TABLE II: Breakdown of systematic errors on partial branch- 
ing fractions B and the ratio R + / . Contributions in the upper 
part of this table are taken to be uncorrelated for B and B + . 





Ae°[io- 2 ] 


AB+[KT 2 ] 


Ai? +/0 [10- 2 ] 


Pmin [ GcV/ c] 


0.6 


0.6 


0.6 


1.0 


Ntags (uncorr.) 


0.126 


0.141 


0.020 


0.020 


Btag background 


0.079 


0.122 


0.014 


0.012 


B -> D 


0.079 


0.041 


0.011 


0.001 


B^D S 


0.050 


0.054 


0.008 


0.002 


X 


0.037 




0.004 


0.006 


D*~ <-» D*° 


0.014 


0.064 


0.004 


0.003 


B — > r 


0.019 


0.020 


0.003 


0.002 


Ntags (con.) 


0.202 


0.253 


0.004 


0.004 


e eff 


0.134 


0.143 


<0.001 


track eff. 


0.085 


0.092 


<0.001 


D,D s ,t^ e 


0.024 


0.020 


<0.001 


conversion, Dalitz 


0.025 


0.039 


0.001 


<0.001 


faked e 


0.020 


0.027 


<0.001 



the semilcptonic decay widths of charged and neutral 
B mesons agree to a precision of 5%, T(B° — ► Xev) / 
T(B+ -» Xev) = 0.998 ± 0.047. 
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